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Abstract 

In this letter we analyze the consequences, for the LHC, of gauge and third family Yukawa 
coupling unification with a particular set of boundary conditions defined at the GUT scale, which 
we characterize as effective "mirage" mediation. We perform a global % 2 analysis including the 
observables M\y, Mz, Gf, a~^, a s (Mz), M t ,m,b(mb), M T , BR(B — > X S ^),BR{B S — > and 
Mfr. The fit is performed in the MSSM in terms of 10 GUT scale parameters, while tan/3 and \i 
are fixed at the weak scale. We find good fits to the low energy data and a SUSY spectrum which 
is dramatically different than previously studied in the context of Yukawa unification. 

PACS numbers: 12.10.Dm, 12.10.Kt, 12.60.Jv 
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Gauge coupling unification in supersymmetric grand unified theories (SUSY GUTs) [1- 
6] provides an experimental hint for low energy SUSY. However, it does not significantly 
constrain the spectrum of supersymmetric particles. On the other hand, it has been observed 
that Yukawa coupling unification for the third generation of quarks and leptons in models, 
such as SO(10) or SU(4) C x SU(2)l x SU(2)#, can place significant constraints on the SUSY 
spectrum in order to fit the top, bottom and tau masses [7-11]. These constraints depend 
on the particular boundary conditions for sparticle masses chosen at the GUT scale (see 
for example, [9, 12-14], which consider different GUT scale boundary conditions). In this 
letter we consider effective "mirage" mediation boundary conditions and show that they 
are consistent with gauge and Yukawa coupling unification with a dramatically different 
low energy SUSY spectrum. The GUT scale boundary conditions are given by an effective 
"mirage" pattern with gaugino masses defined in terms of two parameters, M1/2 an overall 
mass scale and a the ratio of the anomaly mediation to gravity mediation contribution [15- 
17]. Scalar masses are given in terms of (for squarks and sleptons) and mio (for Higgs 
doublets). In addition, the H u and Hd masses are split, either with "Just-So" splitting or 
with a U(l) D-term which affects all scalar masses. Note, as in Ref. [17], we allow for several 
origins of SUSY breaking. For example, the dilaton and conformal compensator fields break 
SUSY at a scale of order M1/2, while the dominant contribution to SUSY breaking is at 
a scale of order m 3 / 2 > mw ~ mw- We fit the low energy observables, Mw, Mz, Gf, ct~^ 
a s (Mz), M t , mt,(mb), M T , BR(B — » X s j), BR(B S — > fi + fi~) and in terms of 12 arbitrary 
parameters. The low energy sparticle spectrum is imminently amenable to testing at the 
LHC. Two benchmark points are contained in Table III. 

Fermion masses and quark mixing angles are manifestly hierarchical. The simplest way 
to describe this hierarchy is with Yukawa matrices which are also hierarchical. Moreover the 
most natural way to obtain the hierarchy is in terms of effective higher dimension operators 
of the form 

45 

W D A 16 3 10 I63+ 16 3 10 — 16 2 + (1) 

This version of SO(10) models has the nice features that it only requires small representations 
of SO(10), has many predictions and can, in principle, find an UV completion in string 
theory. The only renormalizable term in W is A I63 10 I63 which gives Yukawa coupling 
unification 

A = At = Xb = A T = \ Vr (2) 
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at Mqut- Note, one cannot predict the top mass due to large SUSY threshold corrections 

to the bottom and tau masses, as shown in [18-20]. These corrections are of the form 

a 3 M s tan/3 \i A t tan/3 
omb/mb oc % 1 5 h log corrections. (3) 

So instead we use Yukawa unification to predict the soft SUSY breaking masses. In order 
to fit the data, we need 

5m b /m b ~ -2%. (4) 

We take \x < 0, > 0. For a short list of references on this subject, see [7-11, 21-24]. 

We assume the following GUT scale boundary conditions, namely a universal squark 
and slepton mass parameter, m 16 , universal cubic scalar parameter, A , "mirage" mediation 
gaugino masses, 

(where M1/2 and a are free parameters and hi = (33/5, 1, —3) for i = 1,2,3). We consider 
two different cases for non-universal Higgs masses [NUHM] with "just so" Higgs splitting 



m 2 rr = m? n 

H u{d) J-U 



(+)2D (6) 



or, D-term Higgs splitting, where, in addition, squark and slepton masses are given by 

m 2 a = m 2 w + Q a D, {Q a = +1, {Q, u, e}; -3, {L, d}} (7) 

with the U(l) D-term, D, and SU(5) invariant charges, Q a . Note, we take fi, M1/2 < 0. 
Thus for a > 4 we have M3 > 0, Mi, M2 < (i.e. the case of D-term splitting is similar to 
the analysis of Ref. [13]. However our low energy SUSY spectrum is much different.). 

We perform a global x 2 analysis varying the parameters in Table I used to calculate the 
X 2 function in terms of all the observables given in II defined at the electroweak scale as 
discussed in Ref. [14]. We minimize the x 2 function using the Minuit package maintained by 
CERN [25]. Note that Minuit is not guaranteed to find the global minimum, but will in most 
cases converge on a local one. For that reason, we iterate (9(100) times the minimization 
procedure for each set of input parameters, and in each step we take a different initial guess 
for the minimum (required by Minuit) so that we have a fair chance of finding the true 
minimum. We realize that the system is under-constrained and thus we obtain values of 
^ < 1. Nevertheless we believe that the sparticle spectrum is significantly constrained by 
the fit. 
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Sector 


Third Family Analysis 


gauge 

SUSY (GUT scale) 
textures 

SUSY (EW scale) 


1 1 , , . Ala, i ; 
m,i 6 , Mi/2, a , Aq, mio, D 
A 

tan/3, /i 


Total # 


12 



TABLE I: The model is defined by three gauge parameters, aa, Mq, £3; one large Yukawa coupling, A; 6 
SUSY parameters defined at the GUT scale, ttiiq (universal scalar mass for squarks and slcptons), Mi/ 2 
(universal gaugino mass), a (the ratio of anomaly mediation to gravity mediation contribution to gaugino 
masses), mio, (universal Higgs mass), Aq (universal trilincar scalar coupling) and D which fixes the magni- 
tude of Higgs splitting in the case of "Just-so" Higgs splitting or the magnitude of all scalar splitting in the 
case of D-term splitting. The parameters fi, tan j3 are obtained at the weak scale by consistent electroweak 
symmetry breaking. 

Consider first the SUSY spectrum in our analysis. Two benchmark points are given 
in Table III with fixed m 16 = 5 TeV. The first and second family squarks and sleptons 
have mass of order m 16 , while stops, sbottoms and staus are all a factor of about 2 lighter. 
In addition, gluinos are always lighter than the third family squarks and sleptons, and 
the lightest charginos and neutralinos are even lighter. Note the states Xi an d Xi are 
approximately degenerate. In Table III we include the running masses for the chargino and 
neutralino and the dominant one-loop contribution to the mass splitting, AM [30]. Thus 
the chargino signature at the LHC is dominated by the decay x + ~ > X° + 7r+ [31]- This 
typically results in a disappearing charged track since the pion would carry too little energy. 
The present limits from ATLAS are not very constraining [32]. Our LSP is a wino-like 
neutralino. As a result, the thermal abundance of the LSP (obtained using micrOMEGA 
2.4 [33]) is of order 10 -5 due to the large annihilation cross-section to W + W~, i.e. too small 
for dark matter. However, non-thermal production of wino dark matter can give the correct 
abundance [34-37]. Finally, for the two benchmark points, the dominant decay modes for 
the gluino are (calculated using Sdecay [38]) 

• for "Just-so" Higgs splitting - (63% -> x°9 ! 28% ->■ x + bt, -> x~tb and 8% ->■ x°tt) 
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Observable 


Exp. Value 


Ref. 


'- X 3\ 1V1 Z ) 






0.1184 ±0.0007 


[26] 








1/137.035999074(44) 


[261 

L J 








1.16637876(7) x 10~ 5 GeV" 2 


[26] 








80.385 ± 0.015 GeV 


[26] 


M z 






91.1876 ±0.0021 


[26] 


M t 






173.5 ± 1.0 GeV 


[26] 


m b (m b ) 






4.18 ±0.03 GeV 


[26] 


M T 






1776.82 ±0.16 MeV 


[26] 


M h 


125.3 ±0.4 ±0.5 GeV 


[27] 


BR(6 -»• 


S7) 




(343 ± 21 ± 7) x 10~ 6 


[28] 


BR(£ S - 




) 


3.2 x 10" 9 


[29] 



TABLE II: The 11 observables that we fit and their experimental values. Capital letters denote pole 
masses. We take LHCb results into account, but use the average by Ref. [28]. All experimental errors are 
la unless otherwise indicated. Finally, the Z mass is fit precisely via a separate y 2 function solely imposing 
electroweak symmetry breaking. 

• and for D-term splitting - (76% -> x + bi, -> x~tb] 14% -> 3.5% -> x°bb, and the 
rest to light quarks or gluons). 

Note, in the case of "Just-so" Higgs splitting stops are the lightest sfermion, while in the 
case of D-term splitting, sbottoms are lighter. In addition, at low energies, A t , A b are small 
and thus we have small left-right mixing. These affect the gluino decay branching ratios. 
Finally, since both \i and M 2 , M\ are negative and, in principle, we may have light smuons, 
we obtain the correct sign for the SUSY correction to (g — 2) M . However, in practice, our 
sleptons are too heavy to give a significant contribution. In Tables IV and V we give different 
benchmark points, all with x 2 -C 1, in order to present the variation of sparticle masses with 
different values of and M\/2- 

With regards to GUT scale parameters, we find a ~ 12 which corresponds to approxi- 
mately equal dilaton and anomaly mediated contributions to gaugino masses. We also find 
1^3 1 < 1% in the case of D-term splitting or precise gauge coupling unification [39]. 
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In conclusion, we have performed a global x 2 analysis of an SO(10) SUSY GUT with 
gauge coupling unification and top, bottom, r, v T Yukawa unification at Mqut- We have 
analyzed the model for the third family alone. We have shown that the SUSY spectrum 
is predominantly determined by fitting the third family and light Higgs masses and the 
branching ratio BR(B S — >■ fi + 

Our model makes several significant predictions. 

• The first and second family of squarks and sleptons obtain mass of order m 16 , while 
the third family scalars are naturally about a factor of 2 lighter. Gluinos and the 
lightest chargino and neutralino are always lighter than the third family squarks and 
sleptons. 

• Due to Yukawa unification of the third family at the GUT scale we have tan (3 « 50. 
In order to fit the branching ratio BR(B S — > fi + we find the CP odd Higgs mass, 
m a S> Mz- Hence we are in the decoupling limit and the light Higgs is predicted to 
be Standard Model-like. 

• Our LSP is predominantly wino and thus assuming a thermal calculation of the relic 
abundance, we find fl^o ~ 1CT 5 . 

• Xi an d Xi are approximately degenerate. Thus the chargino signature at the LHC is 
predominantly due to the decay x + — > x° + 7r+ - This typically results in a disappearing 
charged track since the pion would carry too little energy. 

• For the two benchmark points, Table III, the dominant decay modes for the gluino are 
for "just so" Higgs splitting - (63% -> x° 9\ 28% x + bi, -> x'tb and 8% 

and D-term splitting - (76% -> x + bi, -> x^tb; 14% -)■ x°tt; 3.5% -> and the rest 
to light quarks or gluons). 
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NUHM 


"Just-so" 


D-term 


77116 


5000 


5000 


y/D 


1877 


1242 


mxo 


6097 


5261 


A 


8074 


5593 


A 1 


-615 


-1294 


M 1/2 


-105 


-100 


a 


11.59 


12.00 


M GUT x 10" 16 


4.50 


2.38 


1 /cx-GUT 


25.11 


25.64 


£3 


-0.039 


-0.007 


A 


0.59 


0.56 


tan /3 


49.43 


48.73 


M A 


1558 


1237 


mi 

11 


1975 


2921 


mr 

01 


2049 


2159 




2473 


3601 


ma 


4905 


5081 


mi 

a 


4944 


4467 




4947 


4477 


Xi 


231.98 


219.11 




232.05 


219.11 


AM = • - .\/ x , 


0.519 


0.438 




882 


874 



TABLE III: Benchmark points and SUSY Spectrum. For each case we have j; 2 < 1 The chargino and 
neutralino masses are tree level and the one loop correction to the mass difference is given by AM. All 
masses are in GeV. 
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mi6 




4000 


4000 


10000 


8000 


Vd 




1725 


1511 


5516 


3207 


mio 




5144 


5079 


13036 


10168 


^4o 




7050 


7542 


15789 


14687 


/' 




-259 


-391 


-1364 


-612 


M 1/2 




-100 


-240 


-120 


-260 


a 




12.00 


11.99 


10.88 


11.58 


Mqut x 10" 


-16 


2.69 


2.27 


2.52 


2.55 


1 / Ct-GUT 




25.29 


25.53 


25.88 


25.76 


£3 




-0.019 


-0.017 


-0.005 


-0.018 


A 




0.616 


0.616 


0.560 


0.606 


tan/3 




50.25 


49.96 


48.68 


49.93 


M A 




1658 


1041 


6975 


2825 


mr 




1308 


1679 


4028 


2751 


mi 




1279 


1760 


3068 


2861 






1613 


1580 


5021 


3282 






3929 


4144 


9659 


7910 


m,- 




3974 


4155 


9876 


7978 


rrie 




3952 


3995 


9808 


7924 






187 


367 


278 


525 






190 


371 


278 


526 


AM 




3.61 


4.54 


0.452 


1.67 


M 5 




858 


1834 


853 


1902 



TABLE IV: Generic Features of the "Just-so" Higgs splitting with the mirage pattern for gaugino 
masses and with different values of miQ and For each case we have x 2 ^ 1- The chargino 

and neutralino masses are tree level and the one loop correction to the mass difference is given by 
AM. All masses are in GeV. 
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mi6 


4000 


4000 


8000 


8000 


\f~D 


1037 


1018 


2531 


1641 


mio 


4598 


4594 


8094 


7351 


^4o 


5588 


5654 


8325 


4810 


A* 


-541 


-591 


-2945 


-2636 


Mi/2 


-100 


-280 


-100 


-280 


a 


12.00 


11.91 


12.00 


10.39 


Mqut x 10" 16 


2.41 


1.87 


1.93 


2.35 


1 / CtQJJT 


25.46 


25.73 


26.05 


26.00 


£3 


-0.011 


-0.009 


0.007 


-0.009 


A 


0.582 


0.599 


0.540 


0.569 


tan /3 


49.20 


49.31 


48.13 


48.70 


M A 


969 


728 


3719 


726 


mr 


2026 


2421 


5178 


5344 


t>l 


1255 


1825 


2634 


4702 




2622 


2644 


5266 


6218 


TTlu 


4091 


4324 


8233 


8105 


mi 

a 


3553 


3825 


6594 


7445 


mg 


3546 


3615 


6607 


7445 


m^o 


215 


529 


226 


529 




216 


531 


226 


529 


AM 


0.554 


2.25 


0.436 


0.475 


M~ g 


867 


2085 


855 


1842 



TABLE V: Generic Features of D-term Higgs splitting with the mirage pattern for gaugino masses 
and with different values of m\§ and M\n- For each case we have x 2 ^ 1- The chargino and 
neutralino masses are tree level and the one loop correction to the mass difference is given by AM. 
All masses are in GeV. 
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